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Cooperativity  at  the  monomolecular  level  binding  of  B  or  C  to  the  common  partner  A  affects 
binding  of  the  other 
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AG'-\  +  AG'\  =  AG°Z+  AGz4 
AAG  =  AGP,-AG\  =  AG:i-AG\ 

C  ase  1 :  AG  ,  >  AG  °4  cooperative  binding 

Case  2:  AG  t  <  AG anticcoperative  binding 
Case  3:  AG\  -  AG:^  independent  binding 


Figure  1  Thermodynamic  cycles  and  cooperativity  (a)  Hypothetical  set  of  bimolecnlar  complexes 
between  component  A  and  two  other  components  ( B  and  C),  with  the  rate  constants,  equilibrium 
constants  and  free  energies  for  complex  formation,  (b)  A  thermodynamic  cycle  for  formation  of  the 
ternary  complex  ABC  by  two  different  possible  routes:  either  B  binds  first,  or  C  binds  first.  There  are 
four  equilibrium  constants  that  describe  the  formation  of  the  various  complexes.  Because  they  converge 
on  the  common  product  ABC,  the  thermodynamics  must  be  independent  of  the  pathway  chosen  around 
the  cycle,  and  constraints  are  placed  on  the  relative  values  of  the  equilibrium  constants  and  hence  the 
free  energies.  The  thermodynamic  coupling  free  energy  (AiC)  gives  the  difference  between  binding  of 
one  component  in  the  presence  of  the  other,  (c)  Definition  of  cooperativity  in  terms  of  binding  of  6  in 
the  presence  or  absence  of  C.  The  two  vertical  binding  reactions  are  gray  to  emphasize  the  comparison 
of  AC"!  and  A<T'4.  It  B  binds  better  in  the  presence  of  C,  the  binding  is  cooperative.  If  B  binds  worse  in 
the  presence  of  C,  the  binding  is  anticooperative.  In  the  third  case,  binding  or  B  is  independent  of  C, 
and  there  is  no  cooperativity. 

Cooperativity  in  macromolecular  assembly 

J  ames  R  Williamson 

volume  4  number  8  August  2008  nature  chemical  biology 


Cooperativity/heterogeneity 


1.  complex  interactions  among  identical 
ligands  binding  to  multiple  sites  on  an 
oligomeric  protein-oxygen  binding  to 
hemoglobin. 

Homotropic  allosteric  regulators— e.g.  02 
Heterotropic  allosteric  regulators— e.g.  2,3 
BPG 


2.  the  thermodynamics  of  macromolecular 
conformational  transitions— protein  folding  or 
nucleic  acid  helix-coil  transitions. 


3.  the  thermodynamics  of  forming  multicomponent 
complexes— multimeric  complexes,  surface  interactions, 
cellular  communication,  organism  organization, 
multicellular  dynamics,  social  structures 

Cooperativity  and  biological  complexity 

Adrian  Whitty  nature  chemical  biology  volume  4  number  8  august  2008 
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Interaction  of  Hemoglobin  with  Three  Uganda:  Organic  Phosphates 
and  the  Bohr  Effect 

(Haldane  cotlli ci«nt /S^-dlphd9>phaxly ct^te/JlnktRe  equatbn^all^tei-io  *(T«t/oiy(feji  binding) 

RUTH  E.  BENESCH  AND  HARVEY  RUBIN 

Department  of  Biocbem  istiy-  Columbia  UiuvttuUi  CoHcro  of  Pliyaifliana  &  Saneom,  New  York,  N .  Y, 
Conuwoiicafc^  tty  Harden  M.  MtCcmndl,  April  16,  lit? 6 


ABSTRACT  The  assumption  that  the  Bohr  coefficient 
(Alog  pse'lM  ApH)  is  equal  to  the  Haldane  coefficient  ( AH  +)  of 
hemoglobin  is  shown  to  be  incorrect  ID  the  presence  of  alien 
Steric  effectors  Such  as  £,3-d Iphos&hogfycerat#,  The  theoreti¬ 
cal  relation  between  th*  two  coefficients  in  the  presence  of 
2,3-diphosphugily  cerate  is  derived.  Experimental  data  on  the 
variation  of  both  coefficients  with  diphosphoglycer* t«  coo- 
central  ion  are  presented  and  shown  to  be  in  agreement  with 
prediction. 


Therefore,  the  liberation  of  diphosphoglycerate  on  oxygen¬ 
ation  must  lead  to  an  increase  in  its  activity  with  increasing 
oxygenation,  and  Eqs.  I  and  2  can  no  longer  apply. 

It  is  the  purpose  of  this  paper  to  examine  the  relation  be- 
tween  the  Bohr  and  Haldane  effects  in  the  presence  of  2,3- 
diphosphoglycerate. 

RESULTS  AND  DISCUSSION 


Box  1  A  timeline  showing  evolution  of  allostery  as  a  concept 

The  timeline  (Fig.  1)  includes  some  of  the  key  exper i merits  and  realizations  of  the  field 
and  some  insight  into  how  the  mentality  of  the  field  has  shifted. 

1903 — The  Bohr  effect  (sigmoidal  binding  curve  of  hemoglobin  to  O2  was  observed). 
1910 — A.  Hill  formulates  the  Hitl  equation  to  describe  the  sigmoidal  binding  of  O2 
to  hemoglobin. 

1958—Rrst  X-ray  structure  (sperm  whale  myoglobin)  solved  by  M  Perutz  and  S11 
J,  Cowctery  Kendrew76. 

1950s— Repression  of  gene  expression,  covalent  modification  of  enzyme  activity,  and 
feedback  inhibition  of  enzymes  are  discovered79. 

1963— J.  Nlonod  renames  regulatory  sues  allosteric  sites ' 

1965 —  J.  Monod,  J.  Ufrman  and  J.-R  Changeux  propose  a  theoretical  model  of  con¬ 
certed  allosteric  transit  ions  (MWC  model)1. 

1966 —  0.  Koshland,  G.  Nemethyand  R.  Films*  propose  the  sequential  model  for  allos¬ 
teric  transitions  tKNF  model)®0. 

1984 — Allosteric  regulation  in  the  absence  of  conformational  change  is  proposed12. 

1 980s— Protein  folding  studies  lead  to  the  concept  that  proteins  exist  in  different 
conformations  ir>  an  "energy  landscape'  81. 

1990s — Mutations,  covalent  mod  if  icatiom  and  changes  in  conditions  such  as  pH  are 
included  as  allosteric  effectors. 

1999 — Allosteric  networks  m  the  PDZ  domain  proposed  by  R.  Ranganathan^. 

2006— Negative  allostehcity  reported  in  the  absence  of  conformational  changes*4. 
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Figure  I  Timeline  showing  I  be  discovery  ami  progression  of  the  concept  of  allosteric  regulation 

in  proteins. 


Logarithmic  scale  of  kcat  and  knon  values  for  representative  reactions  at  25  °C.  The 
length  of  each  vertical  bar  represents  the  rate  enhancement  by  each  enzyme 


ADC  )  arginine  decarboxylase; 

ODC  )  orotidine  5'-phosphate 
decarboxylase; 

STN  )  staphylococcal  nuclease; 

GLU  )  sweet  potato  a-amylase;  FUM  ) 
fumarase; 

MAN  )  mandelate  racemase; 

PEP  )  carboxypeptidase  B; 

CDA  )  E.  coli  cytidine  deaminase; 

KSI )  ketosteroid  isomerase; 

CMU  )  chorismate  mutase; 

CAN  )  carbonic  anhydrase. 


The  Depth  of  Chemical  Time  and  the  Power  of  Enzymes  as  Catalysts 

R.  WOLFE N DEN  AND  M.J  .  SNIDER 
Acc.  Chem.  Res.  2001,  34,  938-945 


How  does  "Biology" cope? 


“After  total  war 
can  come  total  living" 


Mutually  Assured  Destruction:  Cold  War  exhibit  at  the  Smithsonian 


Stringent  response  and  growth  control 


•  Triggered  by  adverse  conditions,  e.g.  starvation 


Transcription  control  (p)ppGpp: 

-  Lack  of  nutrients 

-  Stalled  ribosomes 

-  ppGpp  synthesis 

-  Reprogramming  of  transcription 


Translation  shutdown 

-  Proteases 

-  (p)ppGpp  involved 

-  Activation  of  toxin-antitoxin  modules 

-  Toxin  reversibly  disables  ribosomes 
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Stringent  response  and  growth  control 


•  Triggered  by  adverse  conditions,  e.g.  starvation 


Transcription  control  (p)ppGpp: 

-  Lack  of  nutrients 

-  Stalled  ribosomes 

-  ppGpp  synthesis 

-  Reprogramming  of  transcription 


Translation  shutdown 

-  Proteases 

-  (p)ppGpp  involved 

-  Activation  of  toxin-antitoxin  modules 

-  Toxin  reversibly  disables  ribosomes 


The  Stringent  Response  is  mediated  by  two 
opposing  RelMtb  activities  which  must  be  tightly 

regulated 

1)  pppGpp  synthesis: 

p-p-p-G  +  p-p-p-A<=op-p-p-G-p-p  +  p-A 
GTP  +  ATP  G5  +  AMP 


2)  pppGpp  hydrolysis: 

p-p-p-G-p-p  < - >  PPi  +  p-p-p-G 


pppGpp  alters  RNAP  kinetics  and  mediates  the 
transcriptional  response  to  environmental  conditions  to 
which  Mtb  is  exposed 


The  RAC  Allosterically  Activates  Transferase  Activity 


^ATP  ^GTP  *cat  ^ca/^ATP  ^cat'^GTP 


(mM) 
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RAC  =  RelMtb  Activating  Complex 
Ribosome*Uncharged  tRNA*mRNA 


Heterogeneity  even  within  a  single  molecule 
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Hydrolysis  +  Synthesis 


Hydrolysis  +  Synthesis 
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Figure  I :  Summary  of  RelMth  truncated.  Full-lenght  KclMlb  protein  is  at  the  lop 
followed  by  the  different  truncated  proteins.  Animo  acid  numbers  are  at  the  beginning  and  end 
of  each  fragment  and  corresponding  activity  is  listed  below.  87-187  overlapping  site  is  noted  in 
llie  full-length  RdN1|h 


Cooperativity,  hetergeneity,  stochasticity 
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Another  example:  Controllers  for  nanomachines 
Aerobic  and  anaerobic  respiratory  chain  in  Mtb 


(Rv2194-Rv2196) 


(Rvl451,  Rv2200c,  Rv3043c,  Rv2193 ) 


ATP  synthase 


Aerobic  Pathway 


Anaerobic  Pathway 


Electrons  enter  the  chain  through  NADH  oxidoreductose 

Plot  of  the  NADH-Q2  reductase  reaction  with  varying  Q  concentrations  and  fixed 
concentrations  of  NADH. 

Lineweaver-Burk  plot  [inset),  slopes  (Vmax/Km)  of  the  lines  are  not  affected  by  NADH 
concentration— ping  pong  mechanism 


A 


B 


KmNADH  =  42  uM,  KmQ2=  12.5uM,  Vmax=26  unit  mg-1 


Ping  Pong  tetra-uni  mechanism 

•NDH-2  catalyzes  the  following  two  electron  transfer  reactions: 

•Ndh(Flox)  +  NADH  Ndh(Flred)  +NAD+  (eql) 

•N d h ( F I red )  +  Q  Ndh(FI0x)  +QH2  (eq2) 


A  P  B  R 
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E  is  Mtb  NDH-2,  A  is  NADH  and  B  is  the  quinone 


Phenothiazine  inhibition  of  Mtb  respiration. 


M.  tuberculosis  membranes 

i 

NADU 


2  min 

◄ — ► 


(A)  TPZ  inhibition  of  NADH-dependent 
oxygen  consumption  by  Mtb  membranes 
measured  with  a  Clark-type  oxygen 
electrode.  Respiration  was  initiated  by 
the  addition  of  10  mM  NADH  and  arrested 
upon  the  addition  of  ImMTPZ.  Addition  of 
lOmM  ascorbate  and  1  mM  TMPD 
produced  an  immediate  resumption  of 
respiration. 


Membrane  interaction  domain 
Q338H,  Y364H,  D366G  etc. 


Flavin-binding  domain 
I20T,  G23S,T110A,etc. 


NADH-binding  domain 
F173S,  A178S,  V275E,  etc. 


A  3D  model  of  E.  coli  Ndh  according  to  Schmid  and  Gerloff 
(2004).  Putative  flavin-,  NADH-,  and  membrane-binding 
domains  are  shown  in  ovals. 


A  drug  for  dormant  TB 


Drug 


Rifampin 

Trifluoperazine 

Chorpromazine 

Isoniazid 

Ethionamide 

Capreomycin  sulfate 

Amikacin  sulfate 

Thiacetazone 

Ethambutol 

Streptomycin  sulfate 

p-aminosalicylic  acid 

Ofloxacin 

Tetracycline 

Cycloserine 

Erythromycin 

Dapsone 
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MBC99s  of  17  Drugs  for  Log-phase  and 

6-week-starved  M.tuberculosisH37Rv 
by  cfu  counts. 


We  shall  Destroy  the  respiratory  chain  by 

recognizing  the  heterogeneity  of  its  parts! 
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Type  I:  nuoABCDEFGHIJKLMN 
(Rv3145-Rv3158) 

Type  II:  ndh  (Rvl854c) 

Type  II:  ndhA  (Rv0392c) 


(Rv0534c,  Rv0548c,  Rv0553 ,  Rv0555 , 
Rv0542c,  Rv3853,  Rv0558) 


(Rvl 620 c-1 623c) 


Cytochrome 
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(Rvl451,  Rv2200c,  Rv3043c,  Rv2193 ) 


Aerobic  Pathway 
Anaerobic  Pathway 


phenothiazines 


diarylquinolines 


We  shall  go  on  to  the  end,  we  shall  fight  in  F  ranee,  we  shall  fight  on  the  seas 
and  oceans,  we  shall  fight  with  growing  confidence  and  growing  strength  in  the 
air,  we  shall  defend  our  Island,  whatever  the  cost  may  be,  we  shall  fight  on  the 
beaches,  we  shall  fight  on  the  landing  grounds,  we  shall  fight  in  the  fields  and 
in  the  streets,  we  shall  fight  in  the  hills;  we  shall  never  surrender. 

WSC  J  une  4, 1940 


Can  molecular  computing  say  anything 


based  on  irreversible  nature  of  computation 

The  Fundamental  Physical  Limits  of  Computation 

What  constraints  govern  the  physical  process  of  computing?  Is  a 
minimum  amount  of  energy  required,  for  example,  per  logic  step? 

There  seems  to  be  no  minimum.,  but  some  other  questions  are  open 

by  Charles  H.  Bennett  and  Rolf  Landauer 
Scientific  American  253(1  ):48-56  (July,  1985). 


A  Fredkin  Gate:  Logically  reversible  with  no 
energy  limit  on  the  computation 


A. 

A  B  C  ->  AL  BL  Cl 
111  111 

110  10  1  A 

10  1  110 

1  0  0  1  0  0  B 

0  11  011 

0  10  0  10  C 

0  0  1  0  0  1 
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CAB  is  a  piece  of  DNA  that  we  can  synthesize 
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Figure  2 
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Why  reversible? 


Minimal  energy  expense 

Detection  and  correction  of  intrusion 

Error  checking  by  reversing  computation 
to  recreate  inputs 


Bidirectional  debugging 


In  principle  it  can  take  minimal  energy  to  go 

through  a  biochemical  gate 

DNAn  +  dNTP  DNAn+1  +  PPi 

A  G  =  kt  ln[dNTP/PPi] 

If  dNTPs  are  just  1%  over  the  equilibrium  value: 

A  G  =  kt  ln[10. 1/10]  or  about  0.01  kT 

a  modification  of  an  idea  in  Bennett  and  Landaur's  Sci.  Am 
paper— suggested  using  RNA 


We  synthsized  the  oligonucleotides  and  ran  the  reactions 


Rrsl  Rewiring 
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1st  Step:  Add  I  c*|  AVH 
FG  Output, 

Anneal  &  Rll  in  FQ|1 

C->A  recoders  Cu^t 


2nd  Step 
End  Amplify 

Add  C,  B*  Primers 


3rd  Step: 

Amplify  New  Recoders 
Add  C,  BC  Primers 


1  Y  1 
1st  Step:  Add  I  C*|  Aw|  ET| 
FGa,tt...  FGi2 


AnnealS  Rll  in 
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Output 


2nd  Step: 

End  Amplify 

Add  C,  B’'  Primers 


3rd  Step 

Amplify  New  Recoders 
Add  QBC1  Primers 


Third  Rewiring 
1st  Step:  Add 

FG  Output,  M1M  2nd  Step: 

Anneal  SRII  in  ^  End  Amplify 

G>C  recoders  Output  Add  C,  B*  Primers 


3rd  Step 

Amplify  Rnal  Rewired 
Product 

Add  C,  B  Primers 
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A  fill-in  reaction  occurs  when  the 

0’  value  of  the  recoder  matches 
the  Cw  value  of  the  FG  output. 

This  step  amplifies  properly 
filled  in  recodertoutpwthytrids. 

B7  primers  with  luth  values 
must  be  added  since  the  value 
of  B" is  unknown.  Onlya 
single  C  primer  is  necessary, 
since  tie  value  of  0  is 
determined  by  the  Cf-  A 
recoder  used  in  the  first  step. 

This  value  will  not  change  for 
agven  well  throughout  the 
pocess. 

This  step  generates  new  recoders 
with  a  defined  value  at  A  and 
identical  values  at  B  and  Cw 

The  output  is  a  Cw-B  recoder. 

The  output  of  FGi2  is  added  to 
thenewlymadeC-Brecoders, 
and  tie  rewiring  pocess  is 
repeated. 

This  stepgenerates  C*-C  The  outputof  FG  i2  is  added  to  This  stepamplifies  tie 

recoders  with  a  defined  value  ^  neW|u  ma[je  c^-Crecoders  completely  rewired  input  stand 

at  Aand  B.  and  the  rewiring  pocess  is  '  for  the  next  FG. 

repeated  again. 

Klein,  J  P.,  Leete,  TH.  &  Rubin  H.  A  Biomolecular 
Implementation  of  Logically  Reversible  Computation 
with  Minimal  Energy  Dissipation.  BioSystems  52, 

15-23, 1999. 
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The  gate  works  in  the  lab 
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How  fast  could  one  go  through  one  gate? 

t1/2  annealing: 

3  sec. 

DNA  polymerization  rate: 

15  bases/sec 

For  60  bases  pair  input: 

10  sec 

Some  very  hard  problems  in  nature  (biology- biochemistry)  "solved” 

using  physical  algorithms  that  reduce  the  hardness 

"Problems" 

"Algorithms" 

Search  optimization 

Coopera  tivity 

Hill  climbing— energy  reduction 

Heterogeneity 

Allocation  of  resources 

Stochasticity 

Self  assembly 

Reversible  computation 

Satisfiability 

Controllers  for  nanomachines 

add  your  favorite  problem 
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